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A long-standing problem with the solution of the Schrodinger equation has been its inability to account for 
the electron-positron annihilation in positron hydrogen scattering above the positronium formation threshold. 
This letter shows that this problem has been resolved by the use of the modified Faddeev equations. A number 
of e~e + annihilation cross sections in the energy gap between Ps(ls) and H(n — 2) thresholds are reported 
for both the positron plus hydrogen incoming channel as well as the proton plus positronium incoming channel. 
However the indirect annihilation cross sections after formation of the positronium themselves are well known, 
they will not be included in this report. 

PACS numbers: 36.10.Dr, 34.90.+q 



When the possibility of e + e~ annihilation in flight is ne- 
glected the e + + H collision in the energy region between 
Ps(ls) and H(n = 2) thresholds can have only two open 
channels. They are: the elastic channel and the rearrangement 
channel, where the e + capture the electron from H and form 
a positronium in Ps(ls). Under this assumption annihilation 
can occur only after the formation of the positronium. The 
lifetime of the Ps(ls) is known to depend on the total spin 
S. The singlet para-Ps, 1 Sq, decays into 27 rays with a life- 
time t w 0.13ns. The triplet ortho-Ps, 3 Si, decay into 37 
rays with r w 140ns 0], These well known annihilation will 
not be included in this report. We calculate the e + e~ direct 
annihilation using the standard relation \ 2] 



Van = nro(c/v)Z 
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tq = e 2 /m e c 2 , Z e ff is the annihilation coefficient defined by 
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f is the distance between e + e~ in the three-body wave func- 
tion The calculation of e + e~ in flight thus compensate for 
the shortcoming of the quantum three-body calculation where 
the annihilation channel is closed. 

The modified Faddeev equation yl 01 explicitly separates 
the three-body wave function into the direct channel and the 
rearrangement channel 
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Here x a , y a , a = 1, 2 are the corresponding mass-scaled Ja- 
cobi vectors |4] defined by 



X a = T a (fp - r 7 ) 
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where (aj3j) are cyclic permutations of (123), m a and r a are 
the particle mass and position vectors, 
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FIG. 1: Jacobi vectors in the e + — H system. 



and M = m a + mp + m~ f . 

The Jacobi vectors corresponding to different channels are 
related by the orthogonal transformation 
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%* = (-)^- Q sgn(a-/3)(l-C| Q )V2. 

In this calculation we used the Jacobi vectors given in Fig.^ 
Using bipolar expansion we have 



L=Q M=-L t + X=L 

where a = 1, 2 and the bipolar basis is 
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Y& M (x a ,y a )= J2 (lm\m x \LM)Y^{x a )Y^(ya). 

mi+mx—M 

Here L is the total angular momentum of the three-body sys- 
tem and I, A are the relative angular momenta corresponding 
to the Jacobi vectors x a and y a , respectively. The modified 
Faddeev equations for angular momentum L are given by 
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where 



H, 
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V a is the Coulomb potential, 

w$'A' = (^ M (i Q ,y Q )|v Q |rA' / (^^o)), 

t4 = ^(13) + Vjp(xp,yp) 



V^ s) = V a { a (x a ,y a ) 



V® =V a {l-C a (x a ,y a )), 



C{x,y)) =2{ 1 + exp 
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£0, ?/o> ^ > 2 are the Merkuriev cut-off parameters, and 
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The components tf> a ix(^a, Va) in © are further expanded in 
terms of quintic Hermite polynomial splines for each of the 
variables x a ,y a . The wave functions are solved and normal- 
ized according to the asymptotic wave function |4|,|5[ 

open channels 

^ (a) ~f*+ E K al<J r'. (9) 
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In the Ore gap, the /'s are the product of the standard spher- 
ical Bessel functions and the radial part of the bound-state 
hydrogenic wave functions. K a i a differs from the standard 
if-matrix elements K a i a by only a kinematic factor. 

According to Fig. [0 and Eq. (|2) annihilation takes place 
when X2 = 0. However, it is clear that the integrals in (0 
diverge whenever ip2(x2, IJ2) is involved. Thus, all previous 
calculations of using (0 were limited to energies below 



the positronium formation threshold. Only recently, Ref. |6] 
used an imaginary absorption potential to replace the dynam- 
ics of e _ e + annihilation. In solving the Schrodinger equation 
they found the annihilation cross section below the Ps(ls) 
threshold joins smoothly to the positronium formation cross 
section just above the threshold. But this imaginary potential 
is too week to represent the annihilation dynamics above and 
away from the Ps(\s) threshold. According to this model, 
above the Ps(ls) threshold only indirect annihilation after 
formation is possible. Ref. |7] renormalized the singularity in 
(|2} near the Ps(ls) threshold using the physical 27-lifetime 
of Ps(ls). They also showed the smooth transition of the 
annihilation cross section across Ps(ls) threshold. Theirs is 
essentially a threshold law, valid only near the threshold. Nev- 
ertheless, we will show that the integral 

£inf=<^iK5hM do) 
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is well defined and exists even above the Ps(ls) threshold. In 
principle, the divergent part of (|2ji must be renormalized to the 
physical indirect annihilation cross section after the positron- 
ium is formed. This part is related to the positronium for- 
mation cross sections that we have already reported in 0,01 • 
This report is devoted to the calculation of JlOi . To avoid 
confusion, we define Q0| | to be the annihilation coefficient 
of e~e + annihilation in flight. Thus dlOt truly represents the 



J 



missing open channel, the neglected dynamics in the quantum 
mechanical solution ^ of the e + + H system or other e + + 
atom systems. 

At x 2 = 0, Eq. (jSJl gives 



and 



Vi = -C21/S21X1, y 2 = -I/S2A, 



where Y^f (§2) is the spherical harmonic. Then, 
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where 5 is the phase shift of the scattering problem as a result 
of the normalization according to the asymptotic wave func- 
tion (|9j. 

When the incoming channel is e + + H , tp2 is solely respon- 
sible for all positronium formation cross section. When p+Ps 
is the incoming channel, ipi is responsible for all the hydrogen 
formation cross section . In either case, annihilation in flight 
calculated using J10I exists. Thus Z im f and Z 9m f are calcu- 
lated for these two cases respectively. In Table I we present 
the L = contribution to f° r a number of k values in the 
Ore gap. Table II gives the corresponding annihilation cross 
sections. The point /c = O^la^ 1 is just above the positronium 
formation threshold, where the point k = 0.8612<Iq 1 is near a 
Feshbach resonance. Even though the Estate annihilation is 
relatively smooth in the Ore gap, the sudden increase near a 
Feshbach resonance is quite noticeable. 

All calculations were carried out with a cut-off radii 
y max _ i25a , y max = 228a , except near the Feshbach 
resonance where y max = 150a , y 2 lmx = 250ao was used, 
^inf * s ver y sens iti ve to the cut-off radii y max and y max es- 
pecially when the energies are close to either the lower or up- 
per thresholds. Converged results are obtained only for suffi- 
ciently large y max and y max values. This is demonstrated in 
TablelHIIfor the case k — 0.71, the Ps{\s) formation thresh- 
old is located at k = 0.70653. In contrast, the scattering cross 
sections were well converged at a cut-off radii of less then 
lOOao JUIH], which were cross-checked by another method 
that solves the Faddeev-Merkuriev integral equations by us- 
ing the Coulomb-Sturmian separable expansion approach |9] 
and by numerous previous calculations using the Schrodinger 
equation. 



Calculations with considerable large cut-off radii require 
substantial computer resources, even for relatively low en- 
ergies such as that used in the present calculations where 
all energies are inside the Ore gap between H(n = 2) and 
Ps(n = 1) thresholds. The three-body wave functions were 
obtained by using a quintic spline collocation procedure on the 
Jacobi coordinates |J,|5|]. The dimension of the converged cal- 
culations have reached 175000. Calculations at much higher 
energies and near Feshbach resonances are quite demanding 
in terms of computer resources. Yet such calculations are im- 
portant, for example in the creating of the antihydrogen and 
the studying of its spectroscopy. In the formation of antihy- 
drogen using the process p + Ps — * H + e~, e~e + annihila- 
tion in the outgoing channel can be an important consideration 
when excited states are involved, where numerous Feshbach 
resonance exist. 

Beyond the likely practical applications in antihydrogen re- 
search, this calculation represents the first breakthrough of a 
long-standing problem in positron-atom quantum scattering, 
i.e. its inability to account for the e~e + annihilation above the 
positronium formation threshold. Using Schrodinger equa- 
tion, the calculations in |6] are valid only in the immedi- 
ate vicinity of the positronium formation threshold at k = 
0.70653oq . The modified Faddeev equations enable the ex- 
plicit separation of the e + + H and the p + Ps channels. The 
annihilation cross sections in the e + + H channel are well de- 
fined by JlOb . The singularity involving the p + Ps channels 
are renormalized to the indirect annihilation which can be cal- 
culated using the known physical annihilation cross sections 
of the positronium themselves and the respective positronium 
formation cross sections obtained in the solution of ©, Thus 
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TABLE I: S'-state annihilation parameters in the Ore gap. The wave 
number k's are given in 1 units, where ao is the Bohr radius. 



k 0.71 0.75 0.80 0.85 0.8608 0.8612 
Z x 0.630 0.592 0.531 0.478 0.446 0.403~ 
Z 2 1.699 1.417 0.749 1.055 1.367 7.632 



TABLE II: S'-state annihilation cross sections in the Ore gap. The 
ct's are given in 10 _6 7rao units. 



k 


0.71 0.75 0.80 0.85 0.8608 0.8612 




0.3446 0.3067 0.2578 0.2183 0.2010 0.1819 


02 


0.9293 0.7337 0.3636 0.4820 0.6167 3.442 



a quantum mechanical theory which includes the annihilation 
channel is obtained. 
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TABLE III: 5-state annihilation parameters at k — 0.71 with 6 pairs 
of cut-off radii. 



J/!™* 50 50 55 50 100 125 

yip 3 * 120 150 174 180 228 228 

Zi 0.541 0.574 0.535 0.537 0.633 0.630 

Z 2 0.406 101.8 2.128 0.471 1.766 1.699 



